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Abstract 

The effective potential in the MSSM at the one-loop level is used to evaluate masses of the 
neutral Higgs scalars and to study finite-temperature phase transition. The CP violation 
in the Higgs sector, which is induced by the spontaneous mechanism or by the complex 
parameters in the MSSM through radiative corrections, is determined at zero and finite 
temperatures. 
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1 Introduction 



The scenario of electroweak baryogenesislJ requires that the electroweak phase transi- 
tion (EWPT) to be of first order and that CP violation is effective at that transition 
temperature. The EWPT is of first order only with too light Higgs boson in the mini- 
mal standard model[Q, [|, and it is argued that CP violation in the CKM matrix cannot 
generate sufficient baryon asymmetry. Both the requirements, however, will be fulfilled 
by some extension of the standard model. The minimal supersymmetric standard model 
(MSSM) is one of promising candidates. When some of the scalar partners of the quarks 
and leptons are light enough, the EWPT becomes such a strongly first-order phase transi- 
tion that the sphaleron process decouples just after it with acceptable mass of the lightest 
Higgs scalar Although the masses of the Higgs bosons in the MSSM are constrained 
by some tree-level relations, they receive large radiative corrections from the top quark 
and squarks|5|]. This may broaden a window for successful baryogenesis by the MSSM. 

The MSSM has many sources of CP violations, in addition to the KM phase, such as 
the relative phases of the complex parameters fi, the gaugino mass parameters and scalar 
trilinear couplings, which are effective to generate the baryon asymmetry Q. Besides these 
complex phases, the relative phase of the expectation values of the two Higgs doublets 
could induce the source of baryon number 0. This phase might be induced by radiative 
and finite-temperature effects near the transition temperature || or dynamically generated 
near the bubble wall created at the EWPT, which we call transitional CP violation^. 
This mechanism has been examined dynamically by solving equations of motion for the 
classical Higgs fields connecting the broken and symmetric phases. Then the potential 
for the fields are given by the effective potential at the transition temperature, which are 
approximated by a gauge-invariant polynomial whose coefficients are given by the effective 
parameters at that temperature Jl0|. It was shown that the contributions from charginos 
and neutralinos are important to trigger the CP violation in the intermediate region. 
These analyses contain undetermined parameters such as the transition temperature, 
thickness and velocity of the bubble wall, expectation values of the Higgs fields and the 
magnitude of explicit CP violation at the transition temperature. These quantities should 
be determined by the parameters in the MSSM. 

Now one of our main concerns is whether the EWPT in the MSSM is of first order 
strong enough for the sphaleron process to decouple after it, with acceptable masses 
of Higgs bosons. Recent study of the two-loop resummed effective potential at finite 
temperature suggests that the EWPT is strong enough for 2 < tan /3 < 4, tua > 120GeV 



and a light Higgs with < 85GeVpT[. This and the previous analyses based on the 
one-loop resummed potential did not include the contributions from the charginos and 
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neutralinos. For the parameters which admit the strongly first-order EWPT, we must 
examine whether efficient CP violation exists at the transition temperature. The effective 
potentials used in the previous study of the phase transition were functions of only the 
CP-conserving order parameters, so that they could not evaluate CP violation in the 
Higgs sector at the EWPT. 

In this paper, we use the effective potential in the MSSM, which includes one-loop 
corrections from the top quark, top squarks, gauge bosons, charginos and neutralinos, to 
evaluate the masses of the neutral Higgs bosons and to examine strength of the EWPT. 
The masses are approximated by the eigenvalues of the matrix whose elements are given 
by the second derivatives of the effective potential evaluated at the vacuum at zero tem- 
perature. Although the mass formulas for the neutral Higgs bosons have been found by 
fully-contained one-loop calculations [O and by two-loop calculations [T3J, we adopt this 
method for self-containedness. Without CP violation, the minimum of the effective poten- 
tial is parameterized by by the absolute value and the ratio of the VEVs of the two Higgs 
doublets. Extending the effective potential to include the CP- violating order parameter 
and employing a numerical method to search a minimum of it, we find the magnitude of 
induced CP violation through radiative corrections from the superparticles when some 
of the parameters are complex-valued. To find the transition temperature Tq and the 
magnitude of the VEVs of the Higgs fields at Tq, we numerically calculate the effective 
potential without use of the high-temperature expansion ||14||, and use the same numerical 



method as in zero-temperature case to search a minimum of the effective potential. CP 
violation in the Higgs sector will provide a boundary condition for the equations which 
dynamically determines the profile of the bubble wall [15| . 



This paper is organized as follows. In Section 2, we derive formulas for the neutral 
Higgs boson masses in the MSSM in the absence of CP violation and a numerical method is 
introduced which can be applied to the case with CP violation. The method is extended 
to the effective potential at finite temperature in Section 3. Numerical results on the 
masses of the Higgs scalars and transition temperature and strength of the EWPT are 
summarized in Section 4. Section 5 is dedicated to concluding remarks. The formulas for 
the derivatives of the effective potential are summarized in the Appendix. 

2 Higgs Boson Masses 

If we write the VEVs of the two Higgs doubles as ip^ and <p u , the tree- level potential is 
given by 

V = rn\y\wd + mlvlPu + {[iBip u ip d + h.c) 
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- vivuf + §(^y d )(^u), (2.i) 

where we take rnj = //S to be real and positive by phase convention of the fields. Now 
we parameterize the VEVs as[] 

Then the tree-level potential is expressed as 

Vb = ~™?>i + v\) + im2( V2 2 + v 2 3 ) - m 2 (^ 2 - ^ 4 ) + + «I - < ~ ^lf- ( 2 -3) 

The minimum of this potential is given by 

Vi = v o cos P, v 2 = v sin /3, v 3 = w 4 = 0, (2.4) 

with 

m 2 = m 2 ; tan/3 — -m| cos(2/3), m 2 . = m 2 cot /? + -m| cos(2/3), (2.5) 
so that CP is conserved. The effective potential at the one-loop level is given by 

V cS (v) = V (v) + AV(v), (2.6) 

where 

AV(v) = A g V(v) + A t V(v) + AtV(v) + A x ±V(v) + A x oV(v) (2.7) 
is the sum of the one-loop corrections. Each term is given as follows. 



A g V(v) = 3-2F{m 2 w (v))+3F(m 2 z (v)), (2.8) 
A t V(v) = -4-3-F(m t »), (2.9) 
AtV(v) = 2-3- £ ^HW). (2-10) 



a=l,2 



A x± V(v) = -4 £ ("&(*)). (2- 11 ) 

a=l,2 

A X «V» = -2 E F(mJo(»)), (2.12) 



0=1,2,3,4 



where 



, 9n m 4 / m 2 3\ , 



lr The order parameter Vi can be eliminated by the gauge transformation. So we shall set W4 = except 
when we numerically calculate the Higgs boson masses in the presence of CP violation. 
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which was renormalized by the DR-scheme. The renormalization scale M ren will be taken 
to be the weak scale. The masses-squared in (|2.8| ) and (|2.9|) are given by 



m 2 w {v) 



m 2 {v) 



92 



( 2 , 2 , 2 , 2\ 



m|(«) 



^(v 2 2+ v 2 3 ). 



(2.14) 



m| (f) in fl2.10P are the eigenvalues of the matrix 



M~ 



m 2 ±( 



m 



t L 



+ m 2 t {v) + 3 -^{v 2 - v 2 - v 2 + 
^= [n(vi - iv 4 ) + A t (v 2 + iv 3 )] 



^= \ji(vt + iVi) + A t (v 2 - iv 3 )} 

„2 



m 2 R +m 2 t (v) + g f{v\-vl 



v) in (|2.11 ) are the eigenvalues of the matrix M^ ± M X ± with 



V2 



and m^o(t>) in ( p. 12 ) are the eigenvalues of the matrix M y0 M x o with 



M 2 -^gi{v 2 - iv 3 ) 

g 2 {v 1 - iv A ) -fj, 

t 



(2.15) 



(2.16) 



(2.17) 



M 2 -^2(^1-^4) ^92(v 2 -iv 3 ) ^ 

M 1 ^9i{vi-iVi) -%gi{v 2 -iv 3 ) 

-\92V\ivx -iv 4 ) \gi{vi-iVi) [i 

\ ^92(v 2 -iv 3 ) -^9i(v 2 -iv 3 ) fi 

In general, /i, A t , M 2 and Mi are complex- valued but we assume them to be real until we 
discuss CP violation. 

The minimization conditions of the effective potential relate the mass parameters in 
the Higgs potential to the VEVs of the Higgs fields: 

2 2 a 1 2 , no s 1 9AV(V) 

m, = m,tanp m 7 cos(2p) , 

2 2 a 1 2 , na s 1 9AV(V) 

m 9 = rrin cot a H — m 7 cos[2p) , 

2 3 y 2 z v 11 v 2 dv 2 ' 



;2.18) 



which are used to eliminate m\ and m 2 in favor of vq = |v| and tan/3 = v 2 /v\. These are 
equivalent to the tadpole conditions [0. The masses of the CP-even Higgs bosons are 
the eigenvalues of the matrix 



5gj4ffW d 2 V cff (v) 

dv'f dv\dv2 

d 2 V c ft(v) d 2 V eB (v) I ' 
dvidv2 



[2.19) 



where the derivatives should be evaluated at the vacuum. The mass of the CP-odd scalar 
is given by 

1 d 2 V cS (v) 

(2.20) 



cos 2 j3 dvl 



By use of the second derivatives evaluated at the vacuum are reduced to 

d 2 V eS (v) 2 „ , 2 n d ( 1 dAV(v) . 

i ,2 = ^tan/5 + m|cos 2 /3 + Wl — , 2.21 

d 2 V eS (v) 2 . 2 . 2fl d ( 1 dAV(v) 



m 2 cotf3 + m 2 z sm 2 (3 + v 2 ^ r ( ) , (2.22) 

= -^-^siu^cos^ + ^M, (2.23) 
OV1OV2 OV1OV2 

= m ^ cot ^_i^w + wo. (2 . 24) 

The expressions of the derivatives of AV(«) are summarized in the Appendix. 

In addition to the evaluation of the masses by use of these formulas, we adopt a fully 
numerical method. In this method, the effective potential defined by (|2.6| ) is calculated 
at every v = (^1,^2,^3,^4), where the mass eigenvalues are evaluated numerically. For a 
given set of (vq, tan j3) and m|, the mass parameters in the Higgs potential are determined 
by ( j2.18| ). The minimum of the effective potential is searched by use of the downhill sim- 
plex algorithm [|TB| starting from a randomly generated simplex in the restricted space of 
(v 1, t>2, V3) with t>4 = 0. Once the minimum is found, the second derivatives of the effective 
potential with respect to (vi, v 2 , v$, v±) are numerically evaluated. We have checked, in 
the absence of CP violation, that the minimum coincides with the prescribed (t>o,tan/3) 
and that the four-by-four matrix of the second derivatives is completely divided into to 
the two sectors of CP eigenmodes and the eigenvalues coincides with those obtained by 
use of the formulas above. This numerical method can be applied to the case with CP 
violation. In the presence of CP violation such as the relative phases of complex param- 
eters, CP-violating order parameter V3 is induced and the CP eigenstates of the Higgs 
sector mix to make the mass eigenstates. 

3 Finite- Temperature Effective Potential 



At finite temperatures, the one-loop corrections in ( p.8|) - (p.l2|) are modified to include 
the finite-temperature effects: 



A g V(v;T) = A g V{v) + 



2tt 2 



6/ B (^|+3/ B 



m 2 (v\ 




' m%{v) ' 
2^2 



(3-1) 



A t V(v;T) = A t V(v;T)-12.—I F '-^), (3.2) 



AiV{ V] T) = A { V(v) + 6 ~ X) I* ( ) , (3.3) 
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A X± V>;T) = V^W-^^E^I^I. (3-4) 




A x oV(v;T) = A x oy(«)-2~ S ^l 11 ^)' ( 3 - 5 ) 
where the functions Ib(cl 2 ) and If{cl 2 ) are defined by 

I B:F (a 2 ) = J™dxx 2 log (lTe-^). (3.6) 

The effective potential at finite temperature is calculated at each v by numerically evaluat- 
ing the mass-squared eigenvalues and inserting them into the expressions ( pT5|) . The 
integrals defined in ( |3.6|) are numerically calculated without use of the high-temperature 
expansions] 
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For a given set of parameters, the minimum of the effective potential is searched at 
various temperatures by the method stated in the previous section. Near the transi- 
tion temperature, several numbers of starting simplexes are generated and the minimum 
reached starting from each simplex is found. The temperature at which two degenerate 
minima are found is defined to be the transition temperature Tq of the first-order EWPT. 
Then we examine whether the condition is satisfied for the sphaleron process to decouple 
just after the EWPTfl7|] 

EE = lim MQl > !. (3.7) 

T C T]T C T v 1 

If vq = 0, the EWPT is of second order. We executed this minimum search for various 
sets of parameters to find the order of the EWPT and T c , and measured v c and tan (3 at 
Tq when the EWPT is of first order. 

4 Numerical Results 

4.1 CP-conserving case 

Among many parameters in the MSSM, the mass parameters m\ and m| in the Higgs 
potential are determined by Q2.18Q in the absence of CP violation. Throughout this paper, 



we take v = 246GeV, m w = 80.3GeV, m z = 91.2GeV and m t = 175GeV. The rest of 
the parameters are m|, tan/3, //, A t , m^ L , m t - fl , M 2 and M\. For definiteness, we take 
M 2 = Mi, m h = 400GeV and A t = lOGeV. 

Before presenting the numerical results on the Higgs masses and CP violation, we note 
that the contributions from the charginos and neutralinos are not negligible compared to 
those from the top quarks and squarks. For example, consider the contributions to the first 
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Table 1: Contributions to the equations relating the mass parameters in the Higgs poten- 
tial to the VEVs of the Higgs fields from the stop and charginos. (tan /3, fx) = (2, — 20GeV) 
and (tan/3,/i) = (5, — 50GeV) correspond to the zero stop mixing case. 



tan (3 


At(GeV) 


m t -„(GeV) 


AjmiftGeV*) 


A t ~m 2 (GeV^) 


A x ±mj(GeV 2 ) 


A x ±m 2 (GeV^) 


2 


-20 





-1.323 x 10 3 


-7.402 x 10 3 


5.723 x 10 2 


6.967 x 10 2 


2 


-20 


300 


-1.596 x 10 3 


-1.131 x 10 4 


5.723 x 10 2 


6.967 x 10 2 


2 


-300 





-5.900 x 10 3 


-7.482 x 10 3 


-4.799 x 10 3 


1.263 x 10 3 


2 


-300 


300 


-7.089 x 10 3 


-1.131 x 10 4 


-4.799 x 10 3 


1.263 x 10 3 


5 


-50 





-1.302 x 10 3 


-5.831 x 10 3 


8.050 x 10 1 


7.293 x 10 2 


5 


-50 


300 


-1.573 x 10 3 


-9.017 x 10 3 


8.050 x 10 1 


7.293 x 10 2 


5 


-300 





-4.688 x 10 3 


-5.824 x 10 3 


-3.769 x 10 3 


1.628 x 10 3 


5 


-300 


300 


-5.640 x 10 3 


-8.999 x 10 3 


-3.769 x 10 3 


1.628 x 10 3 



derivatives appearing in (|2.18 ), which have the form of ^^-F' \m 2 ). The contributions 



from the gauge bosons is smaller than those from the top quark, since for tan (3 = 5 

1 dm?(v) 2m? 9 1 dm^Av) q 2 

tK - 1 -y 2 t ^ 1.0526, .p_^ = £1 ~ q.2131, (4.1) 



v 2 dv 2 Vq sin (3 Vi dvi 2 

which multiply — 12F'(mf) and GF'(mw), respectively. For the case of the stop and 



charginos, these factors are replaced with ( | A . 1 3| ) and (|A.22 ) or ( |A.23[ ), which depend not 
only on the couplings but also on /a, tan (3 and the soft-SUSY-breaking masses. When m~ 
is the same order as m x ±, there is no reason for the stop contribution to become much 
larger than the chargino contributions. Now let us denote 

1 dAiV{v) . 2 _ 1 dAiV{v) 



&t m l = a ' A i m 



v\ dv\ " v 2 dv 2 

2 1 dA x ±V(v) 2 1 dA x ±V(v) 

A x ±ml = ir-^A A x ±m 2 = 42 

Vi OVi v 2 ov 2 

and calculate them for M 2 = 300GeV. Several numerical values are presented in Table [l|. 
This suggests that when M 2 ~ ~ m i L i ^t m \ anc i At 777 ! are °f ^ ne same order as 



mf and A^tti^. As seen from ( |A.15|) and ( |A.25| ), the factors multiplying F'(Tn 



2 \ 



which are the couplings squared in the case of the quarks and gauge bosons, are corrected 
by n~ /Rii n xV an d n x^/^x ^ or ^ ne S ^°P S an d charginos, respectively. For the case of 
tan/3 = 5, li = -300GeV and m- tR = in Table [I], nf/FLi ~ 0.53 while n^/R x ~ -13.4 
and w9) /R x ~ 2.85, which are large enough to compensate the difference between the 
gauge and Yukawa coupling constants. As for the neutralino, although its contributions 
cannot be expressed in a compact form as QA.ll ), it is natural to expect the neutralino 



contributions to be the same order as the chargino, as far as M 2 ~ M\. For tan/3 = 5 
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Figure V. Contour plots of the mass of the lighter chargino as function of the MSSM 
parameters M 2 and /i for tan/? = 2, 5 and 20, respectively. All the mass units are GeV. 



and \i = — 300GeV, we have 

= 2.374 x 10 3 GeV 2 , L dA f^ = _ 1M3 x 10 3 GeV 2 , (4.3) 

V\ OV\ t> 2 OV2 

which are the same order as the chargino contributions. As for the second derivatives, 
we find that the contributions from the charginos and neutralinos are the same order as 
those from the stops for the parameters we adopted in the numerical analyses. We also 
observed that if we omit A x ±V(v) and A x oV(v) to determine m\ and m\ in favor of vq 
and tan/? by fl2.18p , the numerical method to search the minimum results in a different 
point in (t> 1; t> 2 )-plane from (v cos /3, v sin 0) with a deviation of about 5% for tan /3 > 10 
and of about 70% for tan/3 < 2. Thus, as long as the soft mass parameters M 2 , Mi and 
mi are of the same order, the contributions from the charginos and neutralinos are 
comparable to those from the top squarks. 

Now we show results on the masses of the neutral Higgs bosons. We examine depen- 
dence of the mass of the lighter scalar m^ on the pseudoscalar mass and M 2 = M\ . In 
practice, we calculate m^ and m& as functions of (m 2 , M 2 ) for a fixed set of (tan f3, /i, mf R ) 
and make a contour plot of m^ in (M 2 , m^)-plane. The mass of the lighter chargino m x ± 
is constrained to be m x ± > 65.7GeV[|11J, which restricts /i and M 2 . According to the 
tree-level mass formula ( |A.21| ), the mass is plotted as a function of M 2 and /i in Fig. [1|. 
This shows that the lower limit is satisfied for the whole range of M 2 we studied, if we take 

> lOOGeV. The limits on the masses of the lighter Higgs scalar and pseudoscalar are 
now mh > 62.5GeV and > 62.5GeV|18|], although more stringent bounds are reported 
m h 75 GeV [15 1 . The results on are plotted in Figs. |^ and [5] for tan/3 = 2, in Figs. [| 
and |5| and in Figs. |B| and [7|. We also calculated m h for = lOOGeV and 200GeV and 
found that behaves in the same manner as the case of |/x| = 300GeV but its value 
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100 200 300 400 500 600 100 200 300 400 500 600 100 200 300 400 500 600 
M 2 M 2 M 2 



Figure 2: Contour plots of as function of the MSSM parameters M 2 and for 
tan/5 = 2, n = -300GeV and m~ tR = 0, lOOGeV and 200GeV, respectively. All the mass 
units are GeV. 



m 




100 200 300 400 500 600 100 200 300 400 500 600 100 200 300 400 500 600 
M 2 M 2 M 2 

Figure 3: The same as Fig. || but for tan/5 = 2, /i = 300GeV and mi R = 0, lOOGeV and 
200GeV, respectively. All the mass units are GeV. 

becomes a bit smaller for smaller |/i|. For /x = 300GeV, in the blank region at small 
vtia and large region, the point v = (vq cos /3, vq sin /3) is not a local minimum but a 
saddle point. This region is broader for larger tan/?, which corresponds to smaller top 
Yukawa coupling. This is because for larger M2, the contributions from the charginos 
and neutralinos to the effective potential, which are negative, dominate over the bosonic 
contributions and make the vacuum unstable. For tan (3 = 2, the experimental bound on 
m,h is satisfied for uia > 200 — 300GeV, depending on \i and M 2 . For tan/3 > 5, it is 
satisfied for ttia > lOOGeV. 

At finite temperatures, the minimum of the effective potential differs from that at zero 
temperature. What we concern here are the order of the EWPT and transition temper- 
ature Tc, at which two minima of the effective potential degenerate, and the location of 
the minimum at Tc when it is of first order. These are important ingredients for the 
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Figure 4: The same as Fig. [2] but for tan/5 = 5, \i = — 300GeV and raj = 0, lOOGeV and 
200GeV, respectively. All the mass units are GeV. 




Figure 5: The same as Fig. || but for tan/5 = 5, /i = 300GeV and m^ R = 0, lOOGeV and 
200GeV, respectively. All the mass units are GeV. 



m t R =0 



m^ R =ioo 



mt R =200 




100 200 300 400 500 600 



100 200 300 400 500 600 
M„ 



100 200 300 400 500 600 



Figure 6: The same as Fig. || but for tan/3 = 20, ji = — 300GeV and m^ R = 0, lOOGeV 
and 200GeV, respectively. All the mass units are GeV. 
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Figure 7: The same as Fig. || but for tan f3 — 20, fj, — 300GeV and mi R = 0, lOOGeV and 
200GeV, respectively. All the mass units are GeV. 



electroweak baryogenesis. For definiteness, we take \x = — 300GeV, M 2 = Mi = 350GeV, 
m~ tL = 400GeV and A t = lOGeV. By use of the numerical method explained in the previ- 
ous section, we search the minimum of the effective potential at various temperatures to 
find the transition temperature. This analysis is done for tan f3 — 2, 5 and 20 with various 
> and two values of m\ being tuned so that ~ 62.5GeV and rrih — 80GeV, 



respectively, for m^ R = 0.0 

For tan (3 = 2 and small rn^. R , it is difficult to have larger than 70GeV as seen from 
Fig. |2|. We have adopted mj = 2.5 x 10 4 GeV 2 , for which ~ 62.5GeV when mi = 0. 
Dependences of vc/Tc, tan /3(Tc) and the masses on are plotted in Fig. §. The 
condition for the sphaleron decoupling after the EWPT ( |3.7| ) is satisfied for m t - B < 75GeV, 
for which rrth < 64GeV and itla — 239GeV. The transition temperature varies from 
T c = 77.3GeV {m- tR = 0) to 88.7GeV (m fjJ = 120GeV) monotonously, tan/3 at T c is 
almost independent of m^ R and remains to be the zero temperature value. 

For tan/5 = 5, we have taken m\ = 3050GeV 2 and 4624GeV 2 , which correspond to 
nih — 62.5GeV and 80GeV, respectively. Tc is monotonously decreasing with respect 
to m- tR and for the former case, 93.2GeV < T c < 100.5GeV, while 93.0GeV < T c < 
100.2GeV for the latter case. Dependence of T c on m§, so that on m^, appears to be 
weak. vc/T c , tan/?(Tc) and the masses are plotted in Fig. § and Fig. |10[ The condition 
( p.7|) is satisfied for m^ R < 50GeV. tan (3 at T c receives finite-temperature corrections to 
become about 20% larger than the zero-temperature value for the case of the larger m^. 

For tan/3 = 20, we adopted m 2 = 2308GeV 2 and 2440GeV 2 . Dependence of T c on m- tR 

and m\ is similar to the previous examples of tan/3 = 5. Now 97.2GeV < Tc < 104.5GeV 

for both choices of m 2 . vc/Tc, tan /3(Tc) and the masses are plotted in Fig. [□] and 

2 For tan/3 = 2 and the rest of the parameters given above, mn cannot be so large as 70GeV. So the 
results are presented only for ~ 62.5GeV. 
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Figure 8: Dependence of vc/Tq (solid curve) and tan /3(Tc) (dashed curve) on mi R for 
m| = 2.5 x 10 4 GeV 2 . rrih (solid curve), rriA (dashed curve) and (dotted curve) are 
also shown for the same parameter set. 
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Figure 9: The same as Fig. [8] but for tan (3 = 5 and m| = 3050GeV 2 
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Figure 10: The same as Fig. |]but for tan/5 = 5 and = 4624GeV 2 
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Figure 11: The same as Fig. || but for tan (3 = 20 and m\ = 2308GeV 2 . 
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Figure 12: The same as Fig. [8] but for tan (3 = 20 and m 2 = 2440GeV 2 



Fig. |12[ In this case, tan/3(Tc) drastically deviate from the zero-temperature value. 

In order to determine the profile of the bubble wall created at the first-order EWPT, 
we must know the global structure of the effective potential at Tc- As an example, we 
show the contour plot of the effective potential at T = and T = Tq for the case of 
tan/3 = 5, m\ 



4624GeV and mf = in Fig. [13]. This shows that two degenerate 



minima at Tc is connected by a valley with almost constant tan (3{Tc) — 6. This is 
also the case for the other sets of the parameters we studied. Some of the baryogenesis 
scenarios based on the MSSM requires that tan (3 varies spatially around the bubble wall. 
But this result implies that tan (3 remains to be almost constant around the wall. 



4.2 CP- violating case 

In the MSSM, there are many sources of CP violation other than the phase in the CKM 
matrix. Among them are relative phases of the complex parameters /i, A t , M2 and M\. 
These also induce CP violation in the Higgs sector, which is the relative phase 9 of the 
VEVs of the two Higgs doublets. This 9 together with the other CP- violating phases 
affect such observables as the electric dipole moment of the neutron. Hence knowledge 
of 9 is necessary to find bounds on the phases of the complex parameters in the MSSM. 
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Figure 13: Contour plots of the effective potential at T = and T = Tc for tan/3 = 5, 
m\ = 4624GeV 2 and m~ tR = 0. 



9 = Arg(f 2 + iva) in the gauge with t> 4 = is determined by minimizing the effective 
potential at T = 0. Even when all the parameters are real, the effective potential could 
have CP-violating vacuum. This is know as the spontaneous CP violation, but in the 
MSSM it inevitably accompanies a too light scalar p0f. We found that this is the case. 
Indeed if we take tan/3 = 5 and mx = and tune m§ = 2113. 3GeV 2 , the effective 
potential has two degenerate minima which correspond to 9 = ±0.318 and tan/3 = 4.917. 
Then the lightest scalar mass is 12.9GeV, which differs from any mass from the mass 
formulas because of large mixing of the CP eigenstates. 

Now we study the effect of explicit CP violation in the complex parameters on the 



CP violation in the Higgs sector. As the first example, we take tan (3 = 5, m^ R 



0, m 2 3 



4624GeV 2 , fi = -300 • e w " and the remaining parameters are set to be the values adopted 
in the previous subsection. For nonzero 5 M , 9 have nonzero value and the scalar and 
pseudoscalar mixes to form the mass eigenstates. By the numerical method, the minimum 
of the effective potential were searched and the second derivatives at the minimum was 
evaluated to calculate the masses of the Higgs bosons, for < 8^ < 0.1. Dependence of 



9, tan {3 and masses of two light bosons on 8^ is depicted in Fig. [HI Within this range of 
8/j,, the derivation of the masses from the values at <5 M is negligible. The induced 9 is the 
same order and has the same sign as 8^. By linearly fitting, we find 9 = 0.8265 • 8^. By 
redefining the fields, we find that the physical CP- violating phase in the mass matrices 
( |2.15| ), ( |2.16|) and (|2.17|) is <5 M + 9. Hence 8^ enhances the magnitude of CP violation. 



As a second example, we put M 2 = Mi = 350 • e tS2 GeV and all the rest parameters 
are taken to be real and set to the same values as the previous example. Dependence of 
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Figure 14: Dependence of 9 (solid curve), tan/5 (dashed curve) and Higgs masses on 




Figure 15: Dependence of 9 (solid curve), tan/3 (dashed curve) and Higgs masses on # 2 . 



9, tan/3 and masses of two light bosons on 5^ is shown in Fig. |T5[ The induced CP phase 
is fitted to 9 = 0.8885 • 62, which has the same sign as the original 62- Since the physical 
CP phase in the mass matrix fl2.16|) is 62 + 9, the CP-violating phase is enhanced by the 
radiative corrections. 

In the scenario of the electroweak baryogenesis, the CP violation around the bubble 
wall is a key ingredient and it is determined by solving the equations of motion with 
the effective potential at Tq- Then the VEVs in the two degenerate minima provide the 
boundary conditions to these equations. Although spontaneous CP violation at T ~ Tc 
could occur, it would accompany a too light scalar at zero temperature. Any way, an 
explicit CP violation is necessary to resolve degeneracy in energy of the CP-conjugate 
pair of the bubble walls, otherwise no net baryon asymmetry would survive the EWPT. For 
the same parameters as the zero temperature case, we plot in Fig. Ill dependences of the 
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Figure 16: Dependence of 9 (solid curve), tan/? (dashed curve) at T = 92GeV. 



induced CP violating phase 9 and tan/3 on the phase of the /z-parameter. The behaviors 
of 9 and tan/3 are almost the same as those at zero temperature, but 9 = 0.8862 ■ 5^. 



5 Discussion 



We have studied the masses of the neutral Higgs bosons and the electroweak phase tran- 
sition of the MSSM by use of the one-loop effective potential. For the parameters we 
adopted, the contributions from the charginos and neutralinos are shown not to be negli- 
gible. We have found that the EWPT is of so strongly first order that the sphaleron 
process after it decouples, for mi < 75GeV when tan/3 = 2, = 62.8GeV and 
m A = 239GeV, for m- tR < 50GeV when tan/3 = 5, m h = 62.8GeV and m A = 70GeV, for 
m t R £ 53GeV when tan/3 = 5, m h = 80GeV and m A = 114GeV, for m- tR < 46.7GeV 
when tan/3 = 20, = 62.7GeV and m A = 62.6GeV, and for mz < 46.8GeV when 
tan/3 = 20, m h = 80GeV and m A = 81GeV. These bounds on almost correspond to 
a bound on the lighter stop mass < m t . For the parameters which permit strongly 
first-order EWPT, we have investigated tan /3 at the transition temperature. It receives 
larger temperature-corrections for larger tan (3 at zero temperature, which corresponds to 
smaller Yukawa coupling of the top quark. This suggests importance of finite-temperature 
contributions from the particles other than the top quark and squarks. We also studied 
CP violation in the Higgs sector, which is characterized by the relative phase 9 of the ex- 
pectation values of the two Higgs doublets. As is well known, the spontaneous mechanism 
to generate 9 accompanies a too light scalar. An explicit CP violation in the complex 
parameters induces 9 of the same order and sign as itself, through radiative and finite- 
temperature corrections. This implies that the physical phases in the mass matrices of the 
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chargino, neutralino and stop are enhanced by the complex phases which are originally 
contained in these matrices. Hence one must take this effect into account to find bounds 
on the parameters in the MSSM obtained from such data as the neutron EDM. 

Some of the mechanism of electroweak baryogenesis in the MSSM requires tan/? to 
vary spatially. But at the transition temperature, it stays almost constant at tan /3(Tc)- 
Then viable scenarios of the electroweak baryogenesis rely on spatially varying 9 and/or 
\v\ in the presence explicit CP violation. The spatial dependence of 9 and \v\ around the 
bubble wall created at the EWPT is examined in |1| . The values of these variables in the 
broken phase at Tq obtained here will serve as the boundary conditions to the dynamical 
equations for (9(x), \v(x)\). These functions in the MSSM will be studied elsewhere |[21||. 

In this paper, we extensively used the one-loop effective potential for self-containedness. 
Extension to the two-loop resummed potential would be straightforward. The two-loop 
resummed potential without the contributions from the charginos and neutralinos yields 
strongly first-order EWPT for a wider range of parameters than the corresponding one- 
loop potential [[nj. We expect that if the higher-order effects are taken into account, the 
effective potential including the contributions from all the particles in the MSSM will 
provide strongly first-order EWPT for a broader region in the parameter space than that 
investigated here. 
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A Derivatives of the Effective Potential 

We present the formulas for the first and second derivatives of the effective potential 
evaluated at a CP-conserving vacuum. Since all the parameters are assumed to be real, 
CP is conserved so that v 3 = v 4 = at the vacuum. But we retain i> 3 to derive the mass 
of the CP-odd scalar according to ( |2.20| ). To see contribution from each species, we give 
the derivatives of the correction to the effective potential from each particle. 



18 



A.l gauge bosons 



The contributions to the effective potential is given by 
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Its first derivatives at the 



1 dA g V{v) 



1 dA g V{v) 
v 2 dv 2 
g| m 2 w (v) 
' 2 32tt 2 



log 



m w( v ) 
M 2 

ren 



gj+gf m%{v) 
2 32tt 2 



log 



m z( v ) 
M 2 

ren 



It is straightforward to calculate the second derivatives: 
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A.2 top quark 

Since the contribution from the top quark (|2 
respect to v\ vanishes. 
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A. 3 top squarks 

The stop contribution is given by (|2.10| ), in which the mass eigenvalues are 
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The first derivatives are given by 
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where the first derivatives of the stop mass-squared evaluated at the vacuum are 
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The general expression for the second derivatives is given by 
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The relevant second derivatives to this expression are 
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A. 4 charginos 
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The derivatives of ( |2.11 ) are evaluated in the same manner as the case of the stop. The 
mass-squared eigenvalues are 
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The first derivatives of the corrections to the effective potential is given by a similar 
expression to ( A.ll ), in which the relevant derivatives of the masses-squared are 
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The second derivatives have the same form as ( |A.16| ) except for the overall coefficient. 
The relevant second derivatives are given by 
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A. 5 neutralinos 



Although the neutralino contribution to the effective potential is given by ( [2.1 2| ) in terms 
of the mass eigenvalues, it is difficult to work out its derivatives, since the eigenvalues 
have complicated forms. To avoid this complexity, we return to the original form of the 
neutralino contribution: 
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where J k denotes the integral over the Minkowskian momentum, the trace is taken over the 
index of 4-dimensional internal space and the spinor index, and D x o is the four-component 
Dirac operator defined by 
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Here the mass matrix M x o is defined by (|2.17 ). The first and second derivatives of (|A.30|) 
have the forms of 



dA x oV{v) 

dvi 

d 2 A x oV{v) 
dvidvj 



- / Tr 

2 Jk 



D x o(k; v 



dvi 



— / Tr 

2 Jk 



dD~J(k-v) dD~hk-v) 



. (A.32) 



22 



The integrand of the first derivative evaluated at the vacuum have the following compact 
form: 
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For the special case of M 2 = Mi which is extensively investigated in this paper, we 
have the following formulas for the relevant derivatives to the masses of the neutral Higgs 
bosons. 
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Here various functions arising from the momentum integrals are defined by 
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where aj = a • '/x. 
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